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Abstract 

The far field acoustic radiation from a lifting air- 
foil in a three-dimensional gust is studied. The acous- 
tic pressure is calculated using the Kirchhoff method, 
instead of using the classical acoustic analogy ap- 
proach due to Lighthill. The pressure on the Kirchhoff 
surface is calculated using an existing numerical solu- 
tion of the unsteady flow field. The far field acoustic 
pressure is calculated in terms of these values using 
Kirchhoff’s formula. The method is validated against 
existing semi- analytical results for a flat plate. The 
method is then used to study the problem of an air- 
foil in a harmonic three-dimensional gust, for a wide 
range of Mach numbers. The effect of variation of the 
airfoil thickness and angle of attack on the acoustic 
far field is studied. The changes in the mechanism of 
sound generation and propogation due to the presence 
of steady loading and non-uniform mean flow are also 
studied. 


N omenclature 


® 1 1 ^2 > ^3 

XuX^X z 


Xc t l 

* 0,*0 


r 


Co 

p' 


physical coordinates 

Lagrangian coordinates of the 

mean flow fluid particles 

modified coordinate for defining <f > 2 

mean flow potential and 

stream function 

steady flow circulation 

speed of sound in the medium 

unsteady pressure 


* Professor, Department of Aerospace and Mechanical 
Engineering, Associate Fellow AIAA. 

1 Research Assistant, Department of Aerospace and 
Mechanical Engineering, Student Member AIAA. 

* Research Scientist, CFD Branch, Member AIAA 


&lj *2.^3 

wave number vector 

A 

Lighthill drift function 

u> 

frequency of the gust 


in radians per second 

<t> 

unknown velocity potential 

<h 

unknown component of <f> 

4> 2 

known component of <f> 


Introduction 


The study of noise generated by propellers and 
turbomachinery blades has received considerable at- 
tention in the recent past. However, much of the at- 
tention has focused on studying the propeller problem 
in its full complexity, and little attention has been paid 
to the more basic problem of noise radiated from an 
arbitrary, lifting airfoil in non-uniform mean flow. The 
characteristics of sound radiated from a lifting airfoil 
in a gust is a fundamental aeroacoustics problem. The 
aim of this investigation is to develop a method to 
calculate the acoustic radiation from a lifting airfoil, 
and to study the effects of thickness and loading. The 
sound intensity and directivity can then be directly 
related to the mean flow Mach number, airfoil geome- 
try and angle of attack, and the characteristics of the 
upstream gust. 

The classical approach in aeroacoustics problems 
is to model the acoustic sources as monopoles and mul- 
tipoles, and to estimate their strength either experi- 
mentally or analytically. Thus, acoustic radiation from 
bodies was mainly studied using Lighthill’s acoustic 
analogy. However, estimation of these acoustic source 
terms from aerodynamic results is difficult and prone 
to error. Even after estimating the sources, one has 
the tedious calculation of the integrals to contend with. 
Clearly, this approach of calculating the acoustics of a 
realistic model is a formidable task. 

The acoustic radiation is the unsteady pressure 
which radiates to the far field due to the unsteady 


1 


pressure fluctuations on the airfoil surface. Recently, 
highly accurate solutions to the problem of an arbi- 
trary, lifting airfoil encountering a periodic gust in 
compressible flow were obtained by J. R. Scott and 
H. M. Atassi [1], [2], [3], [4]. Hence, theoretically 
one can directly calculate the acoustics from the un- 
steady solution given by the Scott- Atassi solver. This 
method, which does not require knowledge of the 
acoustic sources, is the one used in this study. 


Mathematical Formulation 


The current work is concerned with the aeroacous- 
tics of arbitrary lifting airfoils encountering a periodic 
gust. The basic aerodynamic theory was developed 
in References [1] - [5], Since this aerodynamic the- 
ory is used to calculate the airfoil acoustics, we repro- 
duce some of the basic derivation below. This section 
presents a brief outline of the formulation of the un- 
steady boundary value problem. For complete details 
the reader is requested to consult References [3] or [4]. 

The equations governing an inviscid, compressible 
flow past an airfoil at non-zero incidence to a stream 
with uniform upstream velocity Uooi are the continu- 
ity, momentum and entropy conservation equations. 


Dp 

(1) 

H Dt 

(2) 

Ds 

(3) 

— = 0 

Dt 


where, = Jy + U • V, t denotes time, and p, {7, and 
p and s denote the fluid density, velocity, pressure and 
entropy respectively. It is assumed that there are no 
shocks in the flow. Also the fluid is assumed to be an 
ideal, non-heat conducting gas with constant specific 
heats. It is also assumed that there are no imposed 
entropy disturbances or incident acoustic waves. 

Assuming small unsteady perturbations in the 
flow quantities enables linearization of the above equa- 
tions. 

u = 0 o + 5 (4) 

P = Po + p' (5) 

P = Po + p' (6) 

Using the formulation due to Atassi and Grzedzinski 
[5], the boundary value problem is defined as follows. 

A two-dimensional airfoil ( See Figure 1) is placed 
in a potential mean flow with three-dimensional, up- 
stream vortical disturbances of the form 

0 (7) 


The amplitude |a| satisfies \a\ 
defined as 

Xi = UooA 

\T>_ 

* 2 = 


Poo Goo 
A3 = £3 


Uooy and X is 
( 8 ) 
(9) 
( 10 ) 


where, ¥0 is the stream function of the mean flow, x 3 is 
the spanwise spatial coordinate, and A is the Lighthill 
“drift” function [10]. 

The velocity field is expressed as 


U(^t) = U 0 (^ + u(S t t) (11) 


where, Uq(x) is the steady mean flow velocity and 
u(z,<) is the unsteady perturbation velocity. 

The unsteady velocity can be expressed as the 
sum of a vortical part and an irrotational part, such 
that 

u(x,t) = u^ R) + (12) 


where the unsteady potential <t> satisfies the inhomoge- 
neous, non-constant coefficient convective wave equa- 
tion. 


Do, 1 Dq 


(37 ^)--V(poW) 


Dt y c 2 0 Dt 


Po 


= -5-v-(p 0 tr") 

Po 


(13) 


and the unsteady pressure,;/, is given by 


p' - 


(14) 


The unsteady pressure is purely dependent on the irro- 
tational part of the flow and the vortical velocity 
does not induce any pressure fluctuations (to second 
order). 

The boundary conditions on <}> are 


• n = 0 airfoil surface 

(15) 


wake 

(16) 

A[V^ • n\ = 0 

wake 

(17) 

V<j> — ► — as 

X\ — ► — 00 . 

(18) 


Equations (16) and (17) specify the continuity of the 
pressure and normal velocity across the wake, respec- 
tively. For a complete discussion of <£, the reader 
should consult Reference [5]. 

The potential functions <f>\ and <f > 2 are introduced, 
where 

4> — <j>i - 4>2 (19) 
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and fa is a known function which is constructed such 
that 

\<f >2 - fa — * 0 as |x] — ► oo (20) 

where |z| is the distance from the airfoil center. 

The spatial coordinates are non-dimensionalized 
by §; the steady flow density po> by p «>; the unsteady 
pressure p', by Poo^oo|a|; and time, f, by ^7— . The 
normalized wave number, jfci, is also called the reduced 
frequency, and is defined as 


where u and Uoo are the dimensional angular fre- 
quency and the free stream velocity, respectively. 
Henceforth, unless otherwise specified, all quantities 
are assumed non-dimensional. 

The non-dimensional unsteady perturbation pres- 
sure is now given by 


* = /?oo*o (26) 

The unknown potential in the frequency domain is fa. 
We also introduce and fa } where 

<f>2 = (27) 

and 

<P2 = foe''* 1 * (28) 

The aerodynamic code solves for the unknown un- 
steady potential, which will also be referred to as 
the solution to the gust problem. 

The pressure can be expressed in terms of these 
transformed potentials as 

P ' = -Po(-iky + (29) 

We introduce p such that, 


p’ = ~<h ) 


( 21 ) 


where 


2 d 

Dt dt 0 d<t> 0 


The independent variables are chosen to be $o and ^o, 
the mean flow potential and stream function, respec- 
tively. 

Theoretically, the unsteady pressure in the entire 
flow field can be calculated using Equation (21) pro- 
vided the potential functions <j > i and fa are known. Of 
the two, only fa is unknown and is solved for; fa is a 
known function of the upstream disturbance. 

The time dependence of the boundary value prob- 
lem comes entirely through the harmonic term 
A simple change in the dependent variable can trans- 
form the problem from the time domain to the fre- 
quency domain, and significantly simplify the problem. 
The frequency domain transformation is 


p' = pc-«i<+«3*a 


(30) 


and we obtain 

P= ~Po(~iki + t7 o 2 ^)(y>i -¥> 2 ) (31) 

Following the transformations given in Equations (23) 
and (28) for <p\ and y?2i and in Equations (25) and (26) 
for the independent variables, we get 

P = -po[-*(*i + Uo 2 K 1 ){ipi - xl> 2 ) 

+tfo 2 ^(V'i - V- 2 )]e-‘'*‘* (32) 


Formulation of the Kirchhoff Method 


fa = ( 22 ) 

By including the 1T3Z3 term in the transforma- 
tion, the harmonic dependence on the spanwise com- 
ponent X3 is also eliminated, since all of the e** 3 * 3 
terms then factor out from each side of the equation. 

Further simplification is obtained by making the 
following change of both dependent and independent 
variables: 



Vi = 

(23) 

where 

K _ kxMl 
K '- 

(24) 

and 

$ = $0 

(25) 


The Scott-Atassi unsteady aerodynamic solver 
calculates the solution on a grid that extends to the 
far field. In theory it should be possible to compute 
the acoustics using the formulas developed in the pre- 
vious section. However, for symmetric, unloaded air- 
foils and loaded airfoils, the Scott-Atassi code does not 
yield accurate far field results in spite of very accurate 
near field pressure values and lift coefficients (accurate 
within 1% in the validation). This is well illustrated in 
Figure 2. These errors are grid independent since dou- 
bling the grid points on the airfoil from 26 to 52 ( for 
a symmetric, unloaded airfoil with a thickness ratio of 
0.03 ) did not yield any improvement. It is important 
to note that the acoustic energy is a very small fraction 
of the total flow energy. Therefore, a computational 
method such as this requires highly accurate values 
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of the flow potential. However, though the unsteady 
solver yields highly accurate aerodynamic results (near 
field solution), the far field solution is not suitable for 
direct calculation of the acoustics. Since the far field 
solution is very sensitive to the non-reflecting bound- 
ary condition imposed on the outermost grid line, it is 
possible that higher order terms need to be included in 
the Sommerfeld radiation condition. Also the effects 
of numerical dispersion result in reduced accuracy in 
the far field. The high accuracy of the Scott- Atassi 
code in the mid field suggests the use of the Kirch- 
hoff method to circumvent the difficulties arising from 
direct calculation of the acoustics. 

In the Kirchhoff method the unsteady flow field is 
divided into an inner region governed by the inhomo- 
geneous, nonlinear flow equation, and an outer region 
governed by the homogeneous convective wave equa- 
tion. The convective wave equation is transformed to 
the Helmholtz equation in the frequency domain. In 
general, using Green’s theorem, we can calculate the 
value of a function which satisfies Helmholtz’s equa- 
tion, at any point in the outer region, provided we 
know the values of the function and its normal deriva- 
tive on the surface bounding the outer region. If all 
the acoustic sources producing the sound field, and 
all the solid boundaries are confined to a finite region 
R y which is the inner region, enclosed by a Kirchhoff 
surface C, then the acoustic pressure at any point out- 
side C, can be calculated by knowing the pressure and 
its normal derivative on C } and the free-space Green’s 
function for the problem (See Figure 3). 

The Kirchhoff surface is located in the Prandtl- 
Glauert ( ) plane. In the following derivation, the 

superscript V represents the quantity associated with 
the observation point, so that 


r* = 


If (t/>i — fa) satisfies the Helmholtz equation in the 
Prandtl-Glauert plane, then 


where 


and 


(V?. + K 2 )(fa -fa) = 0 

2 _ kjMl _ kj_ 

~ PL PL 

v>.-A+A 

9** 2 dv * 2 


(36) 

(37) 

(38) 


Outside the Kirchholf surface, the non-dimensional 
mean flow velocity is uniformly equal to unity. There- 
fore, 


C = 

where 


-i(ki + f/<j 2 /fi) + Uq 2 


- _ .( *1 . , d 

Co ( /% )+ <9* 0 


Co 


(38) 

(39) 


Also, 

C 0 (V% + K 2 )(fa - fa) = (Vj. + K 2 )C 0 (fa -fa) = Q 

(4°) 

since the operators Co and (VjL-f K 2 ) are linear. Using 
Equations (35) and (38), we arrive at 


(V^^ 2 )P = 0 (41) 

We can now calculate p in the far field using Kirch- 
hoff ’s equation. 

pin = Jm^in?) 


Since the mean flow about the Kirchhoff surface 
is assumed to be uniform, the unsteady perturbation 
potential, and the acoustic pressure satisfy the con- 
vective wave equation. Also since the outer region is 
source free the equation is homogeneous. 

The transformed potential (^i — ^ 2 ) satisfies the 
Helmholtz equation, but p, as defined by Equation (32) 
on the Kirchhoff surface does not satisfy the Helmholtz 
equation. However, Equation (32) can be written as 


p = -po[C(fa - fa)]e ,Ki * 

(33) 

where 

C = -i(k 1 + Uo 2 K x ) + U 0 2 -^- o 

(34) 

We define 


p = -?-e iK '* = C(fa-fa) 

(35) 


Po 


-G(r-|f)|£(r)]da (42) 

where p, the complex conjugate of p, is used in 
Equation (42), because an e tkli dependence was as- 
sumed when deriving the Kirchhoff equation. G is the 
free-space Green’s function for the two-dimensional 
Helmholtz equation, and is given by 


G = - l -4 2 \KR), 


R=\r- rl (43) 


and 


dG 

dn 


.K 




(2) 


(**)[( 




R 


% + <- 


R 


% i <«> 


where r* = ($*,^*) is the observation point in the 
far field, and r = ($,¥) is the “source” point on the 
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Kirchhoff surface. The unit outward normal to the 
Kirchhoff surface is n (See Figure 3). 

It is clear from the foregoing discussion that the 
location and shape of the Kirchhoff surface, C, is of 
little importance, provided it encloses all the acoustic 
sources. However, in practice the choice of C is limited 
by two important considerations. Firstly, it must be 
large enough to include the region of nonlinear behav- 
ior. This is ensured by checking that, the deviation 
in the non-dimensional mean flow velocity values from 
unity, is less than 3% . Secondly, C cannot be so large 
as to lose accuracy of the numerical solution in the mid 
field [4], [7]. Within these limits, the method should 
provide consistent results for small variations in C. 

It was found that the far field results were ex- 
tremely sensitive to the location of the Kirchhoff sur- 
face, even within the aforementioned limits. See Fig- 
ure 4 (b). This problem was traced to numerical errors 
incurred in evaluating the j£ term in Equation (42). 
Since this term involves second derivatives of the solu- 
tion, it could not be evaluated with desired accuracy. 
This motivated us to develop a modified Green’s func- 
tion which eliminated the term altogether. 


The Modified Green’s Function Approach 

In this section we define a new Green’s function, 
G , which is constructed such that it vanishes on the 
Kirchhoff surface, thus eliminating the troublesome 

Gj* term in Equation 42. 

Let 

G(f\r * ) = G(f|r* ) + $(?|r* j (45) 

where G(r|r*) is the familiar free-space Green’s func- 
tion and <?(r|r*) is an unknown function to be deter- 
mined. G is constructed such that 

G = 0 on C (46) 

Also we require G to satisfy 

(V 2 + /C 2 )G(f)r*) = —6(r — f*) (47) 

where r is a variable point in the field, and r* is the 
far field point where we are interested in determining 
p (See Figure 5). The Laplacian operator is 

V 2 = d 2 + d 2 ~ yN 


and Equation (47) implies that 

(V 2 + K 2 )g(r\r*) = 0 (48) 

Also Equations (45) and (46) imply that 

< 7 (r>*) = -G(f|r*) = '-H™{KR) on C 

(49) 

where R = \r — r * |. 

Any solution to the Helmholtz equation can be ex- 
pressed as a summation of Hankel functions, which are 
the fundamental solutions to the Helmholtz equation 
in cylindrical coordinates. 

In order to solve for #(rlr 4 ), we restrict our choice 
of C to a circle of radius t*. This is not a serious 
restriction as the Kirchhoff surface corresponds to a 
£ = constant, line in the computational domain. In 
the Prandtl-Glauert ($-^) domain the level curves of 
the (»/-£) coordinate system tend to rays and circles 
in the far field. To ensure that our approximation is 
valid, we limit the eccentricity of the Kirchhoff ellipse 
to 5%. 

Let r be the observation point on the Kirchhoff 
surface, such that |r| = r* and r* be a “source” point 
in the far field. See Figure 5. For the purposes of 
evaluating the Hankel functions (which are dependent 
only on K and the distance R , between the observa- 
tion point, r, and the source point, r*) the point r 
is equivalent to its “reflection” across the r* vector, 
point rV In other words, the problem is symmetric 
about the vector r*, and one may write 


<;(f|r*) = H% =0 anHW(Kr k ) cos[n(0 - 0*)] (50) 

where only Hankel functions of the second kind are 
considered in order to satisfy Equation (49). 

The coefficients a n are easily obtained as 


a o 


1 


x H?\Kr k ) 


f 


G(r\r*)d</> 


(51) 


and 


where 


*Hj, 2) (Kr k ) 


f G(r]r*) cos[n(0 — 0*)] d(j> (52) 

Jo 


— 0 — 0* and \r\ — r* 


Once </(rlr*) is known, we can rewrite Equation 
(42) as 

Un= J c m^\rldcr (53) 


dG 

dn 


dG dg 
dn + dn 


By definition 

(V 2 -I- K 2 )G(r\r*) = — 6(r — r*) 
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where 


(54) 


~ ~4 Hl ( KR M—r—> dn 


ft ’ 8n J 


and 




(55) 


and we can write 


dft< 2) (ftr t ) _ 
3n 


dH^{Kr k ) dr k <9r t dflt , 

ar t an a« t an J ' ' 

This method gave consistent results for small vari- 
ations in the location of the Kirchhoff surface (See Fig- 
ure 4 (a)) , and was used for all subsequent calcula- 
tions. 

Therefore, using the modified Kirchhoff approach, 
one can obtain accurate far field acoustic pressure val- 
ues from the Scott-Atassi mid field solution. For the 
purpose of numerical calculation, the far field obser- 
vation points are evenly distributed on a circle of suf- 
ficiently large radius ( r is taken to be 50 times the 
half-chord ) in the Prandtl-Glauert ( $ — ¥ ) plane. 
The function g(r]r*) is evaluated in terms of the coeffi- 
cients a n , for the observation point. It should be noted 
that this function depends only on the location of the 
observation point and the Kirchhoff surface, and the 
value of K. Then the derivative is evaluated as per 
Equations (54) and (55). Finally, the far field acoustic 
pressure is obtained using Equation (53), where the in- 
tegral is numerically approximated by a summation of 
the integrand over all the points located on the Kirch- 
hoff surface. Using the relations given in Reference 
9, pp. 28*29, we obtain the corresponding location of 
the observation point in the physical plane. Noting 
that the pressure in the far field goes as the radial 
dependence of the acoustic pressure is eliminated by 
multiplying by y/r. 


Numerical Results 

The modified Kirchhoff method was validated for 
a range of frequencies and Mach numbers against 
the semi-analytical method developed by Atassi and 
Dusey which is detailed in Reference 8. In all the 
cases presented in this paper, the upstream gust veloc- 
ity vector, |a| is of unit magnitude. For the transverse 
gust cases this simply means that = 1. See Figure 6 


(a). In the case of the two-dimensional gust we assume 
that &2 = Jbi, and since the upstream disturbance is 
solenoidal, (a k = 0), we get — a* = a 2 = 0.7071. 
See Figure 6 (b). For the three-dimensional gust cases 
(also referred to as oblique gusts), we assume k 2 = &i, 
and we fix ifc 3 = 0.4. Also we assume ** = — and 
G 2 > 0. These conditions, along with the solenoidal 
condition, (a • it = 0), determine the magnitude of < 33 . 
See Figure 6 (c). 

In Figures 7 to 10, we compare polar plots of the 
acoustic pressure in the physical plane, for selected, 
critical test cases, as obtained by the present approach 
and the semi-analytical approach developed in [8] . The 
length of the line joining any point on the plot to the 
origin represents the magnitude of the acoustic pres- 
sure at a far field point, and the angle the line makes 
with the x-axis is the angular location of the far field 
point. Since the pressure in the far field goes as 
where r is the radial location of the far field point, we 
plot thus removing the dependence on r. 

The transverse gust cases are compared in Figures 
7 and 8 for k\ varying from 0. 1 to 5.0 at a Mach number 
of 0.8. The agreement between the two results is found 
to be excellent. The oblique gust cases are compared in 
Figures 9 and 10 for jfc 3 = 0.4 at a Mach number of 0.8. 
It is found that the agreement is good for ki = 1.0, but 
there is a discrepancy in the magnitude for k\ = 0.5. 

The effects of airfoil thickness, mean loading and 
Mach number on the acoustic pressure emitted from 
an airfoil placed in a transverse or oblique gust are now 
considered. The effect of thickness is first examined, 
and then the effect of mean loading through both angle 
of attack and camber. 


Effects of Airfoil Thickness 

In this section we consider Joukowski airfoil ge- 
ometries of varying thickness. First we study the 
transverse gust cases. The acoustic pressure plots for 
the incompressible flow cases (Moo = 0.1) with k\ 
equal to 1.0 and 3.0 are shown in Figures 11 and 12, At 
low Mach numbers and moderate reduced frequencies 
(ibi < 1.0) the airfoil thickness does not have any sig- 
nificant effect on the acoustic pressure and the dipole 
model is an accurate representation of the airfoil in a 
gust. Hence, those cases are not presented here. How- 
ever, at higher frequencies (see Figure 12) it is found 
that the airfoil thickness tends to tilt the directivity of 
the acoustic pressure toward the upstream direction. 

Figures 13 and 14 show the effect of thickness at 
a Mach number of 0.5. Again we note that the effect 
of thickness is to increase the magnitude of the pres- 
sure in the upstream direction (^ < 0 < 7r). This is 
due to the fact that the incoming flow “sees” a finite 
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rounded leading edge as the airfoil thickness increases, 
as opposed to the infinitesimally thin, flat plate lead- 
ing edge. This causes acoustic propogations in the 
upstream direction. At higher Mach numbers the ef- 
fect of thickness is more pronounced. In Figure 14 it 
is observed that the effect of thickness is to increase 
the acoustic pressure. However, the magnitude of the 
unsteady lift values actually decrease with increasing 
thickness. Hence, a simple dipole model would predict 
a much smaller value of the acoustic pressure. The 
increase in acoustic pressure is due to the effects of 
non-uniform mean flow. Firstly, the non-uniform mean 
flow refracts the acoustic waves radiating from the air- 
foil surface. Even for a uniform mean flow, at high 
Mach numbers, it is known that the acoustic waves 
will propogate at different speeds in different direc- 
tions. This is essentially the Doppler effect. How- 
ever, if the mean flow is non-uniform, the propogation 
speed will depend not only on the direction but also on 
the spatial location of the point under consideration. 
Goldstein, Dittmar and Gelder [9] have shown that the 
sound generated by these “quaulrupole” effects could 
become significant at higher Mach numbers. These 
non-uniform flow effects are discussed in detail in the 
next subsection. 

For the Moo =0.8 case, there is a shock on the 
12 percent thick airfoil and there is as yet no provision 
in the Scott-Atassi code to handle shocks. Therefore, 
only the 3 and 6 percent thick airfoil cases are con- 
sidered at this Mach number. In Figures 15 and 16, 
results are shown for reduced frequencies of 1.0 and 
3.0. The results confirm our previous conclusion that 
airfoil thickness tends to increase the acoustic pressure 
in the upstream direction. This effect is observed only 
at reduced frequencies higher than 1.0. It is also noted 
that the refraction effect of non-uniform flow could ei- 
ther increase or decrease the magnitude of the acoustic 
pressure. 

In Figures 17 through 21, the effect of thickness 
for the oblique gust case is studied. The k\ = 0.5 case 
is not shown at_ M <» = 0.5 since there is no sound 
at that frequency. The parameters should satisfy the 
relation 

*■ > *4 

for sound waves to propogate. At M <*> = 0.5 it is found 
that for the oblique gust case, the reduction in the 
magnitude of the acoustic pressure follows the reduc- 
tion in the unsteady lift . At =0.8, the magnitude 
of the acoustic pressure is almost constant, though the 
unsteady lift values increase. In the following subsec- 
tion, an attempt is made to gain some insight into the 
refraction effects encountered in the cases considered 
this far. This concludes the study of thickness effects. 
In the next section the effects of mean airfoil loading 
are examined. 


Refraction Effects of Non-uniform Flow 

In general it is very difficult to quantise the refrac- 
tion effect of the non-uniform mean flow. However, if 
we take our Kirchhoff surface on the airfoil itself, and 
assume the mean flow to be uniform, we can compare 
the pressure plots and gain some insight into the mech- 
anism of sound propogation in non-uniform mean flow. 

Since the Kirchhoff surface is no longer a circle, 
we cannot use the modified Green’s function approach 
here. Hence, we revert to the free-space Green’s func- 
tion approach. Fortunately it is possible to express 
the normal derivative of the unsteady pressure on the 
airfoil surface in terms of the unsteady velocity on the 
airfoil surface, and the known, steady mean flow quan- 
tities. This enables us to use th*e free-space Green’s 
function approach, which was not possible when the 
Kirchhoff surface was located in the far field. 

The unsteady pressure (physical), is given by 
Equation (21). Differentiating both sides with respect 
to ^ we obtain 


(± + u *£.)?± 

( dt 0 dQ’d* 

.on — ^ = -—(?-) 
< 5 * Vo' 


(57) 


Denoting uj and vj to be the components of the irro- 
tational velocity vector, we can write 


dt _ , 1 V 

a* W ' 


(58) 


d<t> _ , 1 . 

d* { u 0 )vi 


(59) 


Substituting the above in Equation (57), we obtain 


+ 7U 9Uodl._1.9l p^dpo m 

+ 2Uo d*d*~ p 0 d* + p% d* ( } 


Noting that 


and 




u 2 —(—) = 
and substituting for u/ we get 


dU 0 . dU 0 _ 1 dp' p' d Po 

Vl + d* po d* + pi d V 


^ ( 61 ) 
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Identifying |^* with we get an expression for the 
normal derivative of the pressure in terms of the un- 
steady irrotaional velocity and the known mean flow 
quantities. 


w 

dn 


, dU 0 dUo, 

^ Vl W +2ui W + 


P 1 dpo 

Pi d * 


(62) 


On the wake, the unsteady potential <f> is discontinuous 
such that 

<t> + - <f>~ = A(^) 

where denotes the unsteady potential at a point 
just above the wake, and <f>~ denotes the unsteady po- 
tential at a point just below the wake surface. The 
potential jump satisfies 


This formula was used to calculate the normal deriva- 
tive of the pressure on the airfoil surface and subse- 
quently the acoustic pressure was calculated using the 
free-space Green’s function approach. 

The acoustic pressure plots obtained assuming 
non-uniform, and uniform mean flow, for an airfoil 
with a thickness ratio of .12 in a transverse gust with 
jbi = 3.0 at a Mach number of 0.5 are shown in Figure 
22. It is seen that both the magnitude and directiv- 
ity of the acoustic pressure are significantly altered if 
non-uniform flow effects are neglected. In the case of 
uniform flow the magnitude is reduced, and this is as 
expected, since a dipole model of the problem would 
predict acoustic pressure proportional to the unsteady 
lift. 

Similar plots are shown in Figures 23 and 24. In 
Figure 24, a 6 percent thick airfoil in a transverse gust 
with fci = 1.0 and Moo = 0.8 is considered. It is 
shown that dipole models can underpredict the peak 
magnitude of the acoustic pressure by as much as 25% 
in situations where refraction effects are important. It 
is thus seen that the refraction effects can become very 
significant at higher Mach numbers. 


Effects of Angle of Attack 

The effects of airfoil angle of attack on the far field 
acoustic pressure are now examined. The sound radi- 
ated from a 12 percent thick Joukowski airfoil placed 
at varying angles of attack is compared with the sound 
radiated from a symmetric, unloaded 12 percent thick 
Joukowski airfoil. 

In Figures 25 to 30 the incompressible, transverse 
gust cases are presented. The thickness ratio of the air- 
foil is chosen to be 0.12, since the most critical cases 
occur at this value. It must be noted that the angle the 
field point makes with the x-axis in the angle of attack 
cases is (9 - a) and not 9 . Since the only symmetry 
in the problem is about the stagnation streamline, the 
results are presented in this format. The reduced fre- 
quency is varied from 0.1 to 3.0. The most striking dif- 
ference between loaded airfoil cases and the unloaded 
ones is that the acoustic pressure no longer vanishes on 
the wake. This is due to the fact that, in the unloaded 
cases the unsteady potential is an odd function of ^o, 
and hence, the pressure is also an odd function of 


= 0 

For a flat plate or symmetric airfoil, <j> is an odd func- 
tion and we get 

0 + = -<T = 2(A^) 

Consequently, on the wake we get p f to be zero for 
these cases. This is not true for the general case of the 
loaded airfoil. 

At low Mach numbers the effect of the angle of 
attack is to reduce the acoustic pressure on the upper 
side (0 < : (9 — a) : < tt), and to increase it on the 
lower side (x <: (0 — a): < 2 x). For the Jbi = 0.1 
case (see Figure 25), it is observed that there is practi- 
cally no change in the acoustic pressure as we go from 
a = 3° to a = 5°. As the angle of attack is increased 
from 0° to the magnitude of the unsteady lift drops 
by 1.14% . However, there is a drop of only 0.5% from 
a = 3° to a = 5°. This could be one of the reasons for 
the acoustic pressure remaining relatively unchanged. 
It is important to note that the lower lobe of the acous- 
tic pressure actually inreases with increasing angle of 
attack, while the unsteady lift decreases . This implies 
that the effects of mean loading can be important even 
at low Mach numbers. In Figure 26 the single of attack 
is varied from 2° to 10° at a Mach number of 0.1 for 
a reduced frequency of 0.1. It is seen that increasing 
mean loading has the effect of reducing the upper lobe, 
so that at a = 10° the upper lobe is almost completely 
eliminated. 

In Figure 27 the Jbi = 0.5, M^ = 0.1 case is 
presented. The elimination of the upper lobe is now 
observed at even lower angles of attack. For values 
of (0.5 < ki < 1.0), some very interesting effects are 
observed. It is seen that at ki = 0.5, the maximum 
acoustic pressure is along the (0 — a) = direction. 
However, for Jb x = 1.0 (see Figure 28), the maximum 
is along the (9 — a) = ir direction. 

It is of interest to determine whether the large 
peak in acoustic pressure observed in Figure 28, along 
the (9 — a) = x direction, is due to the airfoil leading 
edge thickness or due to mean loading. In Figure 29 
the same case is studied for an airfoil thickness ratio of 
0.03. From this plot we conclude that the peak is not 
due to airfoil thickness. Further parametric studies are 
needed to determine the exact cause of this peak. 

The plot in Figure 30 concludes the incompress- 
ible cases. Here the value of k% is 3.0. For the a = 3° 
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case, the usual suppression of the upper lobe is ob- 
served. However, for the a =5° case, there is a dis- 
tinct rotation of the line of symmetry of the acoustic 
directivity. 

In Figures 31 to 34 the Mqq = 0.5 cases are pre- 
sented. It is observed that at higher Mach numbers 
and low frequencies, the effect of loading is not signif- 
icant. See Figure 31 where jfci =0.1. As the reduced 
frequency is increased from .1 to .5, the effect of an- 
gle of attack is to rotate the acoustic directivity in a 
counter-clockwise fashion (See Figure 32). At k\ = 1.0 
(See Figure 33) there is a very strong effect of angle 
of attack. At this frequency, increasing angle of at- 
tack increases the upper lobe of the acoustic pressure 
by as much as 50% for the a = 5° case though the 
magnitude of the unsteady lift actually decreases by 
1% . Also the maximum occurs at (0 — a) < \ rather 
than at {0 — a) > y, which was documented as the 
effect of thickness in the unloaded case. Furthermore, 
increasing angle of attack introduces a double lobe on 
the lower side of the plot (tc < 0 — a) < 2*). 

Figure 34 concludes the transverse gust cases. It 
is seen that at a reduced frequency of 3.0, numerous 
lobes appear, both above and below the airfoil. There 
are two distinct peaks observed at (6 — o) = j and 
[0 — a) = The region below the airfoil seems to be 
more important at this frequency. 

The oblique gust cases are shown in Figures 35 
and 36. For the ki = 1.0 case, the effect of increasing 
angle of attack is to amplify the acoustic pressure in 
the region above the airfoil, and to suppress the acous- 
tics in the region below. On increasing the reduced fre- 
quency to 3.0, numerous lobes appear. The unsteady 
lift values decrease in magnitude, with increasing angle 
of attack for Jfci = 1.0, while for k x = 3.0 they increase 
with increasing angle of attack. 

The effect of mean loading due to angle of attack 
is to always increase the magnitude of the acoustic 
pressure, and to alter the directivity. The change in 
directivity is specific to the gust, freestream and airfoil 
geometry parameters. It is not possible to generalize 
these trends, rather they have to be investigated on a 
case to case basis. 


Conclusion 

A computational method, based on the Scott- 
Atassi unsteady code, has been developed to obtain 
the acoustic pressure radiated from arbitrary lifting 
airfoils. The effects of airfoil thickness, and mean load- 
ing have been studied for a range of Mach numbers and 
reduced frequencies. It is seen that effect of thickness 
alone is not significant at low Mach numbers and mod- 
erate reduced frequencies. For reduced frequencies less 
than 3.0, the dipole model of the symmetric unloaded 


airfoil is fairly accurate in predicting the acoustics. 
However, at higher Mach numbers refraction effects 
of the non-uniform mean flow become significant and 
the dipole model is not adequate. 

In the loaded airfoil cases, it is observed that the 
dipole model is not valid even for incompressible cases 
with moderate values of the reduced frequency. The 
overall effect of mean loading is to increase the magni- 
tude of the acoustic pressure and to alter the directiv- 
ity pattern. Since the dipole model is no longer valid, 
the magnitude of the unsteady lift can no longer be 
correlated with the acoustic pressure magnitude. 
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SYMMETRIC UNLOADED JOUKOWSXJ AIRFOIL IN A TRANSVERSE CL’ ST 
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KIRCHHOFF METHOD 
DIRECT CALCULATION 
FLAT PLATE 



Figure 2. Comparison between far-field acoustic pressure values calculated using 
the Kirchhoff method and the direct calculation method for a three percent thick 
Joukowski airfoil in a transverse gust at Afoo — 0.1 ; ki = 1.0. The semi-analytical 
results for a flat plate encountering the same gust are also shown. 
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Figure 3. Kirchhoff geometry in the Prandtl-Glauert ( ) plane for the free-space 

Green’s function approach. 



(a). 




Figure 5. Kirchhoff geometry in the Prandtl-Glauert ( $ — 'P ) plane for the modified 
Green’s function approach. 


Figure 4. Comparison between the far-field acoustic pressure as calculated for two 
slightly differing locations of the Kirchhoff surface ( r t = 5.67, r 7 = 5.28 ) using 
(a) the modified Green’s function approach, (b) the free-space Green’s function 
approach for a fiat plate airfoil in a transverse gust at M = 0,8 ; k t = 1.0. The 
semi-analvtical results are also shown. 


ll 


Transverse Gust 


(a). 


U„i 

/i’2 = 0, k 3 = 0, a i — 0, Cl 2 — 1 j Q3 — - 0 



Transverse and Longitudinal Gust 


— ai 

(b). 




1 


A’2 = ki, k 3 = 0, a i = —02, \a\ = 1, a 2 > 0 



Figure 6. Gust conditions for (a) a transverse gust, (b) a transverse 
and longitudinal gust, and (c) a three-dimensional gust. 
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ACOUSTIC PRESSURE PLOT 
FLAT PLATE IN A TRANSVERSE GUST 

1.8 r 

i- 

1 1 J~ KIRCHHOFF SEMI-ANALYTI CaL 



Figure 7. Comparison between the numerically computed far-field acoustic pressure 
and semi- analytical results for a flat plate airfoil in a transverse gust at Af<» — 0.i> , 
t, = 0.1, 0.5, 1.0. 


ACOUSTIC PRESSURE PLOT 
FLAT PLATE IN AN OBUQUE CUST 


KIRCHHOFF SEMI -ANALYTICAL 

KL - 0 5 ' 

K1 > 1 0 -• -• -* 



ACOUSTIC PRESSURE PLOT 



Figure 8. Comparison between the numerically computed far-field acoustic pressure 
and semi- analytical results for a flat plate airfoil in a transverse gust at = 0.8 ; 

i ki = 3.0, 5.0. 


ACOUSTIC PRESSURE PLOT 



Figure 9. Comparison between the numerically computed far-field acoustic pressure 
and semi-analytical results for a flat plate airfoil in an oblique gust at = 0.8 : 
it, = 0.5, 1.0. h = 0.4, |o| = 1. “ = *a = *i. S.k = 0, a, > 0. 


Figure 10. Comparison between the numerically computed far-field acoustic pressure 
and semi-analytical results for a flat plate airfoil in an oblique gust at — 0.8 . 
it, = 3.0, 5.0. = 0.4, |5| = 1, “ = -j, *3 = * 1 . O.k = 0, a, > 0. 
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Figure 11. Effect of thickness on the acoustic pressure radiated by a symmetric, Figure 12. Effect of thickness on the acoustic pressure radiated by a symmetric 
unloaded Joukowski airfoil in a transverse gust. =s 0.1, =* 10 . unloaded Joukowski airfoil in a transverse gust. =» 0.1, k x = 3.0 . 


ACOUSTIC PRESSURE PLOT 



Figure 13. Effect of thickness on the acoustic pressure radiated by a symmetric, 
unloaded Joukowski airfoil in a transverse gust. Af* =* 0.5, k x — 1.0 . 



Figure 14. Effect of thickness on the acoustic pressure radiated by a symmetric, 
unloaded Joukowski airfoil in a transverse gust. M — 0.5, k x = 3.0 . 
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ACOUSTIC PRESSURE PLOT 



ABS(P)*COS(d) 


Figure 15. Effect of thickness on the acoustic pressure radiated by a symmetric, Figure 16. Effect of thickness on the acoustic pressure radiated by a symmetric, 
unloaded Joukowski airfoil in a transverse gust. = 0.8, *, = 1.0 . unloaded Joukowski airfoil in a transverse gust. = 0.8, A, = 3.0 . 


ACOUSTIC PRESSURE PLOT 

SYMMETRIC UNLOADED JOUKOWSKI AIRTOIL IN AN OBLIQUE GUST 
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Figure 17. Effect of thickness on the acoustic pressure radiated by a symmet- 
ric, unloaded Joukowski airfoil in an oblique gust. = 0.5, k, = 1.0. 

i fc 3 = 0.4, |a| = !,«•■ k, = k„ a.k = 0, <Jj > 0. 


Figure 18. Effect of thickness on the acoustic pressure radiated by a symmet 
ric, unloaded Joukowski airfoil in an oblique gust. Af„ = 0.5, t, = 3.0 
k, = 0.4, |a| = 1, “ = *3 = *i. “ * = 0, aj > 0. 
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Figure 19. Effect of thickness on the acoustic pressure radiated by a symmet- Figure 20. Effect of thickness on the acoustic pressure radiated by a symmet- 
ric, unloaded Joukowski airfoil in an oblique gust. Af* = 0.8, k x = 0.5. r ic, unloaded Joukowski airfoil in an oblique gust. M <*, = 0.8, Jb, = 1.0. 


fcj =r 0.4, |aj = 1, jj* = -J, k 7 = k u a.k = 0, a 2 


> 0 . 


* 3 = 0.4, jaj = !,» = k 2 = k u a.k = 0, q 3 > 0. 


ACOUSTIC PRESSURE PLOT 


ACOUSTIC PRESSURE PLOT 




Figure 21. Effect of thickness on the acoustic pressure radiated by a symmet- 
ric, unloaded Joukowski airfoil in an oblique gust. A/» = 0.8. k x = 3 0. 
ib 3 = 0.4, jaj = 1 , 91 = -J, k 2 = k u a.k = 0, a 2 > 0. 


Figure 22. Refraction effect of non-uniform flow on the acoustic pressure radiated by 
a symmetric, unloaded Joukowski airfoil in a transverse gust, thickness ratio=.12. 
= 0.5. k x = 3.0 . 
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Figure 23. Refraction effect of non-uniform flow on the acoustic pressure radiated Figure 24. Refraction effect of non-uniform flow on the acoustic pressure radiated by 
by a symmetric, unloaded Joukowski airfoil in an oblique gust, thickness ratio=.12. a symmetric, unloaded Joukowski airfoil in a transverse gust, thickness ratio=.06. 

= 0.5, Jfc, = 3.0. k 3 = 0.4, \a\ = 1, k 2 = k u a.k = 0, a 2 > 0. .U. = 0.8, *, * 1.0 . 


ACOUSTIC PRESSURE PLOT 



ABS(P)*C0S{*-o) 

Figure 25. Effect of angle of attack on the acoustic pressure radiated by a Joukowski 
airfoil in a transverse gust, thickness ratio=.12. = 0.1, k 2 = 0.1 - 
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Figure 26. Effect of angle of attack on the acoustic pressure radiated by a Joukowski 
airfoil in a transverse gust, thickness ratio=.12. M w = 0.1, k x = 0.1 
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ACOUSTIC PRESSURE PLOT 
LOADED JOUKOWSKI AIRFOIL IN A TRANSVERSE GUST {T/C - 0.12 ) 




Figure 27. Effect of angle of attack on the acoustic pressure radiated by a Joukowski 
airfoil in a transverse gust, thickness ratio=.12. M =0.1, ki = 0.5 . 


Figure 28. Effect of angle of attack on the acoustic pressure radiated by a Joukowski 
airfoil in a transverse gust, thickness ratio=.12, Af« = 0.1, k x = 1.0 . 
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Figure 29. Effect of angle of attack on the acoustic pressure radiated by a Joukowski 
aitfoil in a transverse gust, thickness ratio=.03. A/„ = 0.1. = 10 . 


Figure 30. Effect of angle of attack on the acoustic pressure radiated by a Joukowski 
airfoil in a transverse gust, thickness ratio=.12. Vf„ = 0.1, k, = 3.0 . 
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31. Effect of angle of attack on the acoustic pressure radiated by a Joukowski Figure 32. Effect of angle of attack on the acoustic pressure radiated by a Joukowsk 
in a transverse gust, thickness ratio=.l2. Af. = 0.5, k v = 0.1 . airfoil in a transverse gust, thickness ration. 12. = 0.5. *, = 0.5 . 
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Figure 33. Effect of angle of attack on the acoustic pressure radiated by a Joukowski 
airfoil in a transverse gust, thickness ratio- 12. = 0.5, k, = 1.0 . 


Figure 34. Effect of angle of attack on the acoustic pressure radiated by a Joukowski 
airfoil in a transverse gust, thickness ratio=.12. = 0.5, k t = 3.0 . 
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Figure 35. Effect of angle of attack on the acoustic pressure radiated by a 
Joukowski airfoil in an oblique gust, thickness ratio=.12- = 0.5, k t = 1.0. 
fc 3 = 0.4, | a\ sl.Js -J, Jt, = fc|, al = 0, <z 2 > 0. 
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Figure 36. Effect of angle of attack on the acoustic pressure radiated by a 
■Joukowski airfoil in an oblique gust, thickness ratio=.12. A/oo — 
k 3 = 0.4, ja| = 1, ft = *, = *i, a.k-0, a, > 0. 
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